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The influence of a longitudinal magnetic field on the Coulomb drag current created 
in the ballistic transport regime in a quantum well by a ballistic current in a nearby 
parallel quantum well is investigated. We consider the case where the magnetic field is so 
strong that the Larmour radius is smaller than the width of the well. Both in Ohmic and 
non-Ohmic case, sharp oscillations of the drag current as a function of the gate voltage 
or chemical potential are predicted. We also study dependence of the drag current on the 
' voltage V across the driving wire, as well as on the magnetic field B. 

Studying the Coulomb drag one can make conclusions about the electron spectrum and 
and electron-electron interaction in quantum wells. 
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I. INTRODUCTION 



The influence of a magnetic field on the Coulomb drag are investigated in different geometries. The 
Coulomb drag between two two-dimensional (2D) quantum wells in a strong magnetic field perpendicular 
to the planes of the wells and in the presence of disorder has been investigated in Ref. [1]. fn magnetic 
*-J~ ' field perpendicular to the planes the Hall voltage can be induced in the drag quantum well in the direction 
perpendicular to both direction of the magnetic field and of the current in the drive well [2] , [3] . These 
two geometries can be called transverse. 

The purpose of the present paper is to study the influence of an in- well magnetic field B on the Coulomb 
drag current in the course of ballistic (collisionless) electron transport in a quantum well due to a ballistic 
drive current in a parallel quantum well, fn other words, we consider the longitudinal geometry, i.e. the 
case where the magnetic field is parallel to the applied electric field E and to the plane of the well itself. 

We will concern ourselves with the case of a strong magnetic field that makes the motion of the carriers 
along the field one dimensional and alters the density of electron states. Moreover, we restrict ourselves 
with the quantum limit when only the ground Landau oscillator states are occupied by electrons in the 



two quantum wells, so that 

huj B >fi (f.l) 
Here ojb is the cyclotron frequency while fi is the chemical potential. A theory of electronic transport 



through three-dimensional ballistic microwires in longitudinal magnetic fields at low temperatures has 
been developed in Ref. [4]. Our geometry is similar to that considered in Ref. [4]. However, in the present 
paper we consider much simpler situation of a very strong magnetic field satisfying Eq. (l.f). Later on 
we hope to return to a more general case of a weaker magnetic field where several Landau levels may be 
?3 i involved. 

The magnetic field making the motion of the electrons in the transverse direction one dimensional maps 
the problem under consideration onto the Coulomb drag problem in two one-dimensional wires already 
considered by the authors in Ref. [5] in the Fermi liquid approach. Therefore, our final formulae for the 
Coulomb drag current appear to be similar to those obtained in [5]. Physically the magnetic field may 
play the following important role. It will suppress the tunneling of electrons between the quantum wells 
that, if present, would impede observation of the Coulomb drag. 
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The magnetic field may change the electron quasimomentum relaxation time. Scattering of electrons 
by ionized impurities in sufficiently strong magnetic fields may be even weaker than for B = 0. As for the 
relaxation due to the phonon scattering, a strong magnetic field can alter the density of electron states 
and in the quantum limit the relaxation rate may be bigger than for B — 0. We, however, will assume 
the temperature to be so low that the transport remains ballistic even in the presence of magnetic field. 

We consider the case where the magnetic length 

aB= {w (L2) 

is much smaller than the distance W between the quantum wells of the width L x ~ W each 

W/a B > 1. (1.3) 

This inequality establishes lower bound for the values of the magnetic field for a given distance between 
the quantum wells. For instance, for W <~ 80 nm the inequality requires magnetic fields of the order of 
B <~ 1 T, or bigger. 

It is convenient to break our calculations into several parts. In the first part we will give the principal 
equations of our theory based on the Boltzmann treatment of the transport. We will consider a linear 
response in Sec. III. Next we will discuss a non-Ohmic case in Sec. IV. Comparison of our results with 
the ID Coulomb drag results in the longitudinal geometry and 2D Coulomb drag results for B = will 
be given in Summary. 



II. BOLTZMANN EQUATION 



We consider two parallel quantum wells perpendicular to x axis. The eigenf unctions and eigenvalues 
for a one electron problem in a magnetic field along the z axis in the zth quantum well is [we use the 
gauge A = (0, Bx, 0)] 



^OpyP* = n—i- ^ ~ exp(ip y y/fi + iPzz/K), (2.1) 

y/L y L z \ a B J 

«*-"< + ! r + &- (2 ' 2) 

Here m is the effective electron mass, ipo [{% — x Py)/ a B~\ is the wave function of a harmonic oscillator in 
the ground state oscillating about the point x p = —a B p y /fr = —py/mus- The wave function ipop y p z 
describes a state for which the electron probability distribution is large only within the slab of the width 
w as symmetrically situated about the plane x = x Py and falls off exponentially outside the slab. As 
we consider the case L x >• we will assume the wave function to be equal "4>op y p z if x Py is within 
the quantum well and zero otherwise. In what follows we will need the matrix elements of the functions 
exp(± iqr) between two stationary states. We have 

(Qpyje^'lQPwP*) - e ±lqAxp ^'y )/2 
xe-B^/ 4 e- (l "'- !B 't )a/4o2fl <y piJ ..±R„V sJ .»±ft 9 »- (2-3) 
The diagram representing Coulomb drag effect is illustrated in Fig. 1. The external driving force enters 
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FIG. 1. Coulomb drag diagram. Here the labels 2, (1) stand for the drive (drag) quantum wells. 



the diagram through nonequilibrium distribution function represented by the solid lines marked by the 
symbol 2 indicating that they represent the drive quantum well. 

Now we embark on analysis of the conservation laws for the collisions of electrons belonging to two 
different quantum wells. We have 

e (1) +e (2) -e (1) +e (2) (2 4) 

where eg,? = U h2 + Hlj b /2 + p 2 z /2m. 
The solution of Eq.(2.4) is 

^Qz=p' z -Pz- (2.5) 
The 5— function describing energy conservation can be recast into the form 

5(eil + efX - e^ z+hqz - e™_ n J = ^[Uq z - (p' z - „,)]■ (2-6) 

Therefore, the initial quasimomenta p z and p' z after the collision become p z + fiq z — p' z and p' z — hq z =p z , 
i.e. the electrons swap their quasimomenta as a result of collision. 

Following [5,6] we assume that the drag current in the quantum well 1 is much smaller than the drive 
ballistic current in the quantum well 2 and calculate it by solving the Boltzmann equation for the quantum 
well 1. We have 



Vz 



OA F$ y (p z ,z) 
dz 



-/(12){ F (1) )F (2) }! (2J) 



where i^ 1,2 ) are the electron distribution functions in the quantum wells 1 and 2 respectively, and the 
collision integral I^ 12 ^ {F^-\ F^} takes into account the interwell electron-electron scattering 

/( i 2 ) {jF (i) 5F ( 2)}= <;:Sr 2p '-^(^, fe ,^, P ,, P ;)5. (2.8) 

P'zP'yl'^ 

In this expression the sum over p' y should be determined by the requirement that the a:-center of the 
oscillator function is within the second quantum well. The requirement imposes the constraint 
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'w + L.-^)<^<^(w + L. + ^), (2.9) 
and the product of distribution functions <S is 

* - (l ~ (l - F$,_ hq ) - ^ + ^U (l - *£>) (l - . (2.10) 

*g! - - Ms') + 0[-f*]/(eop, - Ms*) + AF« (2.11) 

where 

r , (1 for »2>0 
9[V ' ] = \0 for v z <0- 

Here we assume that the electrons move ballistically within the quantum well and the electrons moving 
from the left and right reservoirs have the chemical potentials fi B L = pb — eVJ-j/2 and fj, B R = ps + eV^j/2 
respectively. We introduce also the drag voltage induced across the drag quantum well due to the 
quasimomentum transfer from the driving quantum well, i.e. we assume an open circuit for the drag 
quantum well. 

The solution of Eq.(2.7) is (here we omit the equilibrium part) 

AFgl=-( Z ±±)±lW{ F V,FW}, for P >>°> (2.12) 
Using the particle conserving property of the scattering integral 

J2 /(12) {F (1) )F (2) } = Q (213) 

VzPy 

we get for the total current in the drag quantum well defined as 

J =fE^ 5 (2-14) 

Z PvVz 

the result 

J=-e ]T I {12) {F {1 \F^} + e±- £ v z [f(e 0pz - - f(s 0pz - Ms fl )]. (2.15) 

Py,(Pz>0) Py,(P^>°) 

In these equations the sum over p y is restricted by the requirement that the cc-center of Landau oscillator 
must be within the quantum well, so that —HL x /2a B < p y < fiL x /2a B . Introducing the density of 
states (including spin) per unit quasimomentum interval 



N(p z )d Pz = 2 tn _l^" J v d Pz (2.16) 

we have 



2eV d HL x L y f°° ( df(e-p B ) 



For the degenerate electron gas this expression can be written as 

e 2 L L 

^Ohm = -^^f^d- (2-18) 
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Here the number of Larmour circles covering the cross section of the quantum well L x L y /2ira B appears 
instead of the number of open channels in the ID situation. 

We assume that only the ground Landau oscillator state is occupied, so that 

1 3 

Ux + -Hujb < Ms < Ui + —hu>B ■ (2.19) 

Taking into account Eq.(2.6) we obtain for the Coulomb scattering probability Eq.(2.8) 

X (0p zPy \e-< x -^y-^ z \0p z + hq zPy + hq y }(0p z + fiq zPy + Tiq y \e^ x+ ^y + ^ z \0p z p y ) 
X(0p> z p' y \e^ x+ ^y + ^ z \0p' z - hq z p' y - hq y )(0p' z - hq z p' y - hq y \e~^ x '^y'^ z \0p' z p' y ) . (2.20) 

Here we use the unscreened Coulomb potential postponing discussion as to when this approximation can 
be justified until the last section. 

To calculate the drag current we iterate the Boltzmann equation in the interwell collision term that 
we assume to be small. Therefore one can choose the distribution functions in the collision term to be 
equilibrium ones, e.g. = /(4p — Ms) f° r the nrst quantum well. 

We assume, in the spirit of the approach developed by Landauer [7], Imry [8] and Biittiker [9] the 
drive quantum well to be connected to reservoirs which we call 'left' I and 'right' r. Each of them is 
in independent equilibrium described by the shifted chemical potentials Ms = Ms — eV/2 and \i r B = 
Mb + eV/2, where Ms is the equilibrium chemical potential in the magnetic field. Therefore, the electrons 
entering quantum well from the 'left' ('right') and having quasimomenta p' z > (p' z < 0) are described by 

-^Op' = f( £< 0p' ~ Mb) [Fop 1 — Z(4p' — Mb)] an< ^ we see that the collision integral Eq.(2.8) is identically 
zero if the initial p' z and final p' z — q quasimomenta in the drive quantum well are of the same sign. This 
means that only the backs cattering processes contribute to the drag current. 

Due to Eq.(2.6) we are left only with p' z < (since we are restricted according to Eq.(2.15) by the 
constraint p' z — Kq z = p z > 0) and obtain in view of the <5-function in Eq.(2.20) the following product of 
distribution functions in the collision term 

-p - *S!*S! r - <) - *£!') - (i - *S!) (i - F Sl r ) . ( 2 - 21 ) 

or 

V = f(e$ x - Mb)/(4S - Mb)[1 - /(4pl - Mb)][1 - /(4pl " Mb)] 

-/(4p1 - Mb)/(4p1 - Mb)[1 - /(4pl - Mb)][1 - /(4pl - Mb)]- (2-22) 
This equation will be analyzed in the following sections. 

III. LINEAR RESPONSE 

In this case eV/T <C 1 (we assume the Boltzmann constant to be equal 1) and Eq.(2.22) can be recast 
into the form 

-P = ^/(4pl - Mb)/(4p1 - Mb)[1 - /(4pl - Mb)][1 - /(4pl - Mb)]- (3.1) 
Shifting the integration variable p' y — > p' y + h (W + L x )/a 2 B we have for the drag current 
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drag 



eV 2ir ( Air e 2 



t n 



1 

2^h 



2L z dp z 



2ttH Jq 2-kTl (p z + p' z ) 



hL x /2a 



B 2L v dp y dp' v 



.900 [(Pz + P'z) IK, (Py -Py)/ h ] 



where 



A(k Z7 k y , q y ) 



hL,j2a 2 B 2lT h 2tt?T 
9oo(k z ,k y ) = J ^-e- a B<A 2 {k z ,k y ,qy), 



2tt 



(3.2) 



(3.3) 



(3.4) 



(kl + ql) 

1± = lx + The last integral over the centers of the Larmour circles (since x p — —a B p y /h) in Eq.(3.2) 
plays the role of an effective Coulomb interaction potential between the electrons freely moving along the 
direction of applied magnetic field. 

hL x /2a 2 E 



I 



2Lydp y dp'y 
2irh 2^h 



9oo 



Pz+p' z Py-Py 



h 



h 



-fiL x /2a| 



2Ly 



L x /a\ 



dk v ky < goo 



Pz+P z h L _ h 

h ' a% \ 



+ 3oo 



P^ + Pz b _h. 
h ' v al 



(3.5) 



We keep only the first term in this expression since the second term includes a faster oscillating exponent 
~ exp (iq x (W + L x )) as compared to the oscillating exponent in the first term <~ exp (iq x W). 
As W/a,B >lwe can sufficiently simplify the expression for 500 • We obtain 



-^-A 2 {k z , k y , q y ). 



9oo(k z ,k y ) = e a z* 

j —00 ^ n 

dq x eilxKW+L^-alky] e -\W+L x -a 2 B k y \y/$+kl 



A(k z ,k y ,q y )^ J ^ _ , fr? 



2 J«2 + kl 



Finally, the interaction term acquires the form 

r hL x /2a 2 B 2Ld dp > 



hL x /2a% 27r h 2-kTi 



9oo [kz,(Py-p' y )/ti] 



where 



For a > 1 



4a B (27ra_Bfc z ) 3 



e a B k *<f>(2Wk z ) 



$(a) 



(3.6) 
(3.7) 

(3.8) 
(3.9) 
(3.10) 



This result for the effective interaction (3.8) can be explained as follows: Larmour circles within the 
quantum wells of the width L x ■ l/(k z L x ) near the surfaces contribute to the interaction. The number 
of interacting circles from two quantum wells is 

L y L x ■ 1/ (k z L x )^ 



a B 



a B 
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The sum over q x ,q' x <~ k z ,q y ~ \Jk z jW leads to a factor (k z L x ) 2 ■ L y \Jk z jW . The exponential decay of 
the drag with the distance between the quantum wells W is a consequence of one-dimensional character 
of the drag in the strong longitudinal magnetic field. Combining all these factors and multiplying the 
result by U 2 ~ (4?r e 2 ) 2 /(L x L y L z ) 2 k* we arrive at Eq.(3.8). 

The product of the distribution functions in Eq. (3.2) is a sharp function of p z and p' z at small tem- 
peratures, acquiring nonzero values only at p z ,p' z ~ Pp ± T/vp. We assume that the quasimomentum 
interval T/vp is much smaller than h/W 

T « M. (3.11) 

Here we wish to note that the Boltzmann treatment of transport phenomena requires that the uncer- 
tainty in longitudinal momentum must be smaller than the same momentum interval h/L z <C T/vp. 
These two requirements automatically lead to the inequality W <C L z . We assume that the last inequality 
holds. 

According to our assumptions we can regard the interaction term in Eq. (3.2) as the slowly varying 
function and obtain 



J drag J ° A£ b £ b\ 2 T 



-p c p 

where 



sinh 



Ul2 

2T 



l -2 



(3.12) 



^(Airhy a B {a B kpY K ' 



Here we introduced notations TJn — U\ — Ui and mvp = Pp = \Jlm\\iB — C/i — Hub/ 2], kp = Pp/Ti. 
We assume that the electrons remain degenerate in the magnetic field 

e | = MB -^-^ » T. (3.14) 

We consider the quantum limit, i.e. the case when all electrons belong to the first Landau level 

ef < huj B . (3.15) 

Since the electron concentration Nb under this condition is related to the chemical potential by the 
equation 

mfiuiBPp 

7T 2 h 3 

Eq.(3.14) and Eq.(3.15) lead to 



Nb = =F (3.16) 



T«f)!<^ = ^. (3.17) 
2m m Tilob 

The first inequality in this relation is weaker than Eq.(3.11) if Sp ~ Uujb and Wkp > 1 . Introducing 
the electron concentration N and the chemical potential fi for B = given by 

N=^l, s F = ,-U 1 (3.18) 



one can rewrite Eq.(3.17) as 



N B \ ( s F \ 
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Note that the second inequality in this expression does not depend on the electron mass and can require 
magnetic fields stronger than Eq.(1.3) [thus imposing a constraint on the electron concentration, or, if 
the latter is given the inequality may require stronger magnetic fields than is required by the Eq.(1.3)]. 
For instance, in a magnetic field of the order of B ~ 10 T the electron concentration N must be smaller 
than 2.7 • 10 17 cm" 3 . 

Considering the case of the aligned quantum wells, so that U\ = U2 [otherwise the effect is exponentially 
small, cf. Eq.(3.12)] and putting N — Nb we obtain 



drag- J o 4T y £p J 



eV (T\ 2 (3hLJB^ 4 
2e F 



(3.20) 



Jo = - 



e 5 mL y L z kp 



2s F 



e 12(2 £F /3^ B ) 3 $ mkp 



2e F 



(3.21) 



9K 2 (477^)3 

The drag current is a rapidly increasing function of the applied magnetic field, as the latter increases the 
density of states and decreases the transferred Fermi momentum. 
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FIG. 2. Drag current versus dimensionless magnetic field b — TxubIef for two values of the interwell distances 
W — 40nm (1) and W = 50nm (2). Other parameters are given in the text. 



To make an estimate of the current we put m = 0.07m e , Hlub ~ £f — 14mcV, k — 13, L z ~ L y 
W = 40nm. 
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In the linear response regime we can introduce a drag resistance, i.e. we can introduce the coefficient 
that depends only on the quantum wells parameters and relates the drive current JA r i ve in the quantum 
well 2 to the induced voltage in the drag quantum well Jd Y [ ve RY) = ^d- Here the drive current in the 
quantum well 2 is 



drive 



= —V 



e 2 L x Ly 
2wh Tra 2 B 



(cf. Eq.(2.18)) and is determined by the condition of zero total current J = ^j ra g 
the drag quantum well in Eq.(2.15). 



Ji 



Ohm 



R D = 



irh Eb T L z 



e 2 e F e F L y (k F L x ) 2 \ Ae F 



e 12(2e F /3^ B ) 3 $ 



4Wk F 



2 e F 



3 Tiloi 



(3.22) 



in 



(3.23) 
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where we introduced effective Bohr energy Eb — me 4 /K 2 h 2 and 

Pf = YW^^ Pf ' PF = V2me F . (3.24) 
With the given above parameters we have the following estimate for the transresistance 

R D ~ 0.4 mfi. 

Now let us discuss when one can neglect the screening. Since the transferred momenta are q z <~ 2pf/h 
we may not take into account the screening of the Coulomb potential if the inverse screening length is 
much smaller than the transferred momentum. We estimate the screening length at a transferred energy 
- T as 



1 



'ire 2 NB £ 



ln^. 
T 



The required inequality can be written as (we put Nb = N) 



N i/3 e 2 



In 



£F I £F 



T \ Uuir 



< £ F 



£F 

TlLOB 



We will assume this inequality to be satisfied. 



IV. NON-OHMIC CASE 



(3.25) 



(3.26) 



The product of distribution functions Eq.(2.22) can be recast into the form 



V = 2 sinh (eV/2T) exp{(e« - fi B )/T} exp{(e$ - ^ B )/T} 



(2) 



(4.1) 



x/(4V - - Ms - eV/2)f(e { p )> - fi B )f{e^ - M b + eV/2) 

As V is a sharp function of p z and p' z one can take out of the integral all the slowly varying functions 
and get 

r HL x /2al 2Lydpy dp/y 



I 

JO 



2 Z {Pz+P' z ) J-hL x /2a 2 B 2lT h 2TTU 

eV_Ui2 

L y m 2 a%T 2 exp (2a B k* f ^wffi sinh ( £^ 4T 2T 



goo [(pz +P z )/h, (p y -p'y)/h] 



4(47r?i) 3 (a B fcf) 6 
The drag current is 

1 



eV 
4T 



U12 
2T 



J. 



drag 



= Jo 



T 



sinh 

\2T 



2T 



eV 
4T 



, . eV U 12 

sinh 

4T 2T 



sinh[^ + ^ 
4T 2T 



(4.2) 



Ul2 

2T 



eV 
4T 



U12 
2T 



sinh 



eV 
4T 



Ul2 

2T 



sinh 



eV 
4T 



Ul2 

2T 



(4.3) 



For eV <C T one gets from Eq. (4.3) the result of Eq. (3.20). Let us consider the opposite case eV ^> T. 
In this case one gets a nonvanishing result for Eq.(4.3) only if | C12 1 < eV/2 and one obtains the following 
equation for the drag current 



'drag 



Jo 



(eV_\ _ (U12 
\A£ F J \2e F 



itlLOB 

2e F 



(4.4) 



Thus the drag current vanishes unless eV > 2|J7i 2 |. 
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V. SUMMARY 



We have developed a theory of the Coulomb drag between two quantum wells in a strong longitudinal 
magnetic field. We have considered a comparatively simple limiting case where only the lowest Landau 
level is occupied. The strong magnetic field makes transverse motion of an electron one dimensional. 
These one dimensional electron states can be visualized as quantum "tubes" or "wires". Therefore, the 
Coulomb drag problem in this situation becomes similar to the Coulomb drag problem between two 
parallel nanowires. 

It is interesting to compare our results with two different geometries of experiment. First, let us 
consider the influence of magnetic field on ID Coulomb drag for the longitudinal geometry. In this case 
the magnetic field is directed along z axis and is parallel to ID nanowires. For simplicity, we assume that 
the confining potential in the absence of the magnetic field is 

U(x,y) = ^(x 2 +y 2 ). (5.1) 
The applied magnetic field shortens the radius of the state qb so that it becomes 



\A + (B/B, 



2mfl 



where do is the radius in the absence of the magnetic field. For the lowest Landau level we have 

11 fi 2 v 2 

■exp(-p 2 /44), £p = - r + ^. (5.3) 



'27T as 2ma„ 2m 



B 

The wave function of the electron in the second wire can be obtained by a gauge transformation of the 
wave function in the first one. Since the interaction term g is not phase sensitive we are left only with a 
shift by the distance W between the centers of the wires in the argument of the wave function (5.3). As 
a result, one gets for the interaction 



9(2p?) = 46-"'*'^. 



n 2 



(5.4) 



where Iq{x), K (x) are the modified Bessel functions. The quasimomentum 



must satisfy the inequality 



P B f = 



h 2 



T<<(p " )2/2m< ^ < 5 " 5 > 

since we have assumed that only the lowest Landau level is occupied. Here Nf is the electron density 
per unit length in magnetic field. The expression (5.4) demonstrates that provided the magnetic field 
goes up the localization radius a# of the wave functions suppresses the probability of the backscattering 
processes. Note that if one assumes = Nl, where Nl is the electron density per unit length for 
B = then the effective interaction depends on the magnetic field only via a b ■ Therefore in this case the 
magnetic field does not change the magnitude of transferred momentum, in contrast with the previous 
case where such a change leads to a rapid increase of the drag current in a strong magnetic field. The 
drag current is 
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Joi = -^^L z k B F g{2p B F ). (5.7) 

Second, we can compare our results with the drag between two two-dimensional quantum wells [10] in 
the field-free case. In this case the transresistance P12 is proportional to 

pi2 ~ T \kldY {klar (5 - 8) 

where k$ is the single-quantum well (two-dimensional) Thomas-Fermi screening wave vector of the order 
of the inverse effective Bohr radius, d is the interwell distance. First we note that the temperature depen- 
dence of the Coulomb drag between two (three-dimensional) quantum wells in the strong magnetic fields 
is weaker than for the drag in two dimensions. Second, we note that in the latter case the contribution 
from the backscattering processes can be neglected as compared to the small angle scattering contribu- 
tion with transferred momenta < q < 1/e? <C k$ while in our case only the backscattering processes 
are important (this is again a consequence of one-dimensionality of the Coulomb drag problem in the 
quantum limit). 
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